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Abstract 

We solve the Kadanoff-Baym equations for nonequilibrium initial configurations 
of the </) 4 -theory in 2+1 dimensions and compare to explicit solutions of general- 
ized transport equations for the same theory. The latter transport equations are 
derived from the Kadanoff-Baym equation in a first order gradient expansion in 
phase space and explicitly retain the off-shell dynamics as inherent in the time- 
dependent spectral functions. The solutions of these equations compare very well 
with the exact solutions of the full Kadanoff-Baym equations with respect to the 
occupation numbers of the individual modes, the spectral evolution as well as the 
chemical equilibration process. Furthermore, the proper equilibrium off-shell distri- 
bution is reached for large times contrary to the quasiparticle Boltzmann limit. We 
additionally present a direct comparison of the solution of the generalized transport 
equations in the Kadanoff-Baym and Botermans-Malfliet form; both solutions are 
found to agree very well with each other. 
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1 Introduction 



Nonequilibrium many-body theory or quantum-field theory has become a major topic of 
research for transport processes in nuclear physics, in cosmological particle physics as well 
as condensed matter physics. In particular the dynamics of heavy-ion collisions at vari- 
ous bombarding energies - as representative of complex systems with strongly interacting 
particles or fields - has always been a major motivation for research on nonequilibrium 
quantum many-body physics. The initial state of such collisions represents an extreme 
nonequilibrium situation, whereas the final state might exhibit a certain degree of ther- 
malization in central reactions. Indeed, at the presently highest energy heavy-ion collider 
experiments at the Relativistic Heavy-Ion Collider (RHIC), where one expects to create 
experimentally an intermediate state of deconfined matter (denoted as quark-gluon plasma 
(QGP)), there are experimental signatures for an early thermalization accompanied with 
the generation of a very large pressure. Such a rapid thermalization cannot be explained 
by binary on-shell parton scattering using cross sections from perturbative QCD (pQCD) 
PP . Presently open questions in high-energy nuclear physics are, if i) the quantum off-shell 
equilibration proceeds much faster than ordinary two-body on-shell Boltzmann kinetics 
and if ii) multi-parton interactions play a decisive role. 

We here primarily focus on the first issue and present solutions for the full quantum 
time evolution of nonequilibrium 4 -theory in 2+1 dimensions in comparison to general- 
ized off-shell transport equations - derived in the last years [21 El as well as the on-shell 
quasiparticle Boltzmann limit. The scalar 4 -theory, of cause, is much simpler than a 
non-abelian gauge-field theory like QCD, however, provides a suitable testing ground for 
approximation schemes to the field theoretical many-body problem. 

We recall that a proper basis for nonequilibrium quantum many-body physics is pro- 
vided by the closed-time-path (CTP) functional integral ('Schwinger-Keldysh' [U El HE 
IB]) formalism. The essential merit of this real-time Green function technique is its 
applicability to all kind of quantum nonequilibrium phenomena. In this context the most 
prominent role is played by the causal Dyson-Schwinger equations of motion for the one- 
particle Green functions, i.e. the Kadanoff-Baym (KB) equations P ITUl ITT| IT^l ITS} IT3]. 
The latter, furthermore, have served as a basis for the derivation of generalized transport 
equations, which play the central role in almost all practical simulations of nonrelativis- 
tic and relativistic nonequilibrium many-body systems. In the past, a major impetus 
has been to derive semiclassical Boltzmann-like transport equations within the standard 
quasiparticle approximation ^3^21 El QH|- More recently, off-shell extensions by means 
of a gradient expansion in the space-time inhomogeneities have been formulated for sev- 
eral systems: for a relativistic electron-photon plasma for transport of electrons in 
a metal with external electrical field [2D] and in semiconductors ^HHH], for transport of 
nucleons in heavy-ion reactions at intermediate energies [22] 5 for transport of particles in 
4 -theory ^21 I2H] as well as for transport of partons in high-energy heavy-ion reactions 

[23 123 123 123- 

Within all these derivations the particle self-energies - entering the KB equations - 
have to fulfill consistency relations in order to preserve symmetries of the fundamental 
Lagrangian |2H1 I2H]- in this respect an elegant way for achieving selfconsistency is 
provided by the two-particle- irreducible (2PI) effective action [T21 EDI EI]- Including 
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the background field and the correlation function on the same footing, collisions can 
be naturally included by taking into account contributions beyond the leading mean-field 
order of the generating functional. Due to the construction in terms of full Green functions 
this method ensures the validity of all essential conservation laws. 

Moreover, in the nuclear physics context the quantum dynamics of the spectral func- 
tion is a lively discussed issue in the microscopic modeling of hadronic resonances with 
a broad mass distribution. This is of particular relevance for simulations of heavy-ion 
reactions where e.g. the A-resonance or the p-meson already show a large decay width 
in vacuum. Especially the p vector meson is a promising hadronic particle for show- 
ing possible in-medium modifications in hot and compressed nuclear matter (see e.g. 

EH1 El]), since the leptonic decay products only interact electromagnetically. Hence, 
a consistent formulation for the transport of extremely short-lived particles beyond the 
standard quasiparticle approximation is needed. On the one side, there exist purely for- 
mal developments starting from a first order gradient expansion of the underlying KB 
equations EI3 ESI E7J EH] , while on the other side already first practical realizations for 
various questions have emerged El EH HOI HI] ■ However, a fully ab initio investigation 
for a strongly interacting field theoretical system in 3+1 dimensions -without any further 
approximations - does not exist so far as well as a rigorous demonstration of the reliability 
of the different schemes for off-shell transport dynamics. 

Although the Kadanoff-Baym equation has been of wide interest, (numerical) investi- 
gations of its full solution are rather scarce for strongly interacting systems. First studies 
in the nuclear physics context have been performed by Danielewicz |42j . who investigated 
the initial distribution of two occupied non-relativistic Fermi spheres in momentum space 
for spatially homogeneous systems. In comparison to a standard on-shell semiclassical 
Boltzmann equation the full quantum Kadanoff-Baym solution showed a larger collective 
relaxation time which was attributed to quantum interference and off-shell effects. Similar 
quantum modifications in the equilibration and momentum relaxation have been found 
in J3] and for a relativistic situation in Ref. [43J. Furthermore, the spectral properties 
of nuclear systems in the nonrelativistic regime have been first examined on the basis of 
the KB equations in Refs. |H[ H5] . We mention that solutions of quantum transport 
equations for semiconductors ^3] US] - to explore relaxation phenomena on short time 
and distance scales - have become also a very active field of research (cf. Ref. [UJ for a 
recent overview). 

The basic aim of our present study is to examine exact solutions of the KB equations 
for strongly interacting (scalar) fields in comparison to the solution of generalized kinetic 
equations as obtained from a gradient expansion of the full Kadanoff-Baym equation. 
Although such generalized transport equations have already been used in the past, a 
check against the full solution is still lacking. To this aim we will investigate the scalar 
4 -theory for the case of 2+1 space-time dimensions since - first of all - it superceeds the 
rather artificial structure of the 4 -theory in 1+1 dimensions, that has been investigated 
before. In the latter case on-shell scattering processes are (col-)linear and thus do not 
lead to an equilibration in accordance with an appropriate change of the momentum 
distribution. Thermalization in 1+1 dimensions is only mediated by off-shell scattering 
processes. This is a rather unusual situation, since it is known from the investigation of 
heavy-ion collisions by means of quantum kinetic theories, that binary on-shell scattering 
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plays an important role. In the 2+1 dimensional case both, on-shell and off-shell collisions 
influence the equilibration behaviour. Thus the theory in 2+1 dimensions is a well-suited 
laboratory in order to investigate the differences between the full Kadanoff-Baym theory 
and approximate schemes as given by the on-shell Boltzmann limit or by generalized 
off-shell transport equations. Furthermore, the divergence structure of the 4 -theory in 
2+1 dimensions allows for a controlled renormalization procedure because only two self- 
energies - containing pure vacuum parts of Green functions - diverge in the ultraviolet 
regime. Thus the renormalization can be performed equivalently to the standard zero 
temperature vacuum case. This is different in 3+1 dimensions, where complete classes of 
self-energies are divergent, which incorporate non- vacuum parts of the Green functions as 
well. In the latter case more involved renormalization schemes should be used analogous 
to those constructed for self-consistent approximations at finite temperature [JSJ HHl ED] • 
Our work is organized as follows: In Section 2 we will present the relevant equations 
for the nonequilibrium dynamics in case of the 4 -theory and provide the Kadanoff-Baym 
equations within the three-loop approximation of the 2PI effective action. Furthermore, 
we briefly review the derivation of generalized transport equations and discuss different 
limits. Section 3 is devoted to numerical studies on equilibration phenomena within the 
full Kadanoff-Baym theory (cf. Ref. [ST]) in comparison to the generalized transport 
equations in either Kadanoff-Baym (KB) form [5] or Botermans-Malfliet (BM) approx- 
imation [22J. We will also address the differences in the solutions of the two different 
versions (KB and BM) of the off-shell transport theory and explore the validity of the 
gradient expansion employed. We close this work in Section IV with a summary of our 
results and an outlook on future extensions. 

2 Derivation of Generalized Transport Equations 

In continuation of our previous study in Ref. [HI] we investigate the out-of-equilibrium 
quantum dynamics of interacting scalar fields by the example of the 4 -theory that allows 
to test theoretical approximations in a variety of aspects jH3 EH1 EH EH] • Its Lagrangian 
density is given by (x = (t,x)) 



where m denotes the 'bare' mass and A is the coupling constant determining the interaction 
strength of the scalar fields. 

For the derivation and investigation of generalized transport equations we first have 
to specify the necessary notations and the explicit form of the Kadanoff-Baym equations. 
Though the details of the derivation are given in Ref. |2] we briefly recall the main steps 
in order to keep the paper self-containt. 

2.1 The Kadanoff-Baym Equations 

A natural starting point for nonequilibrium quantum theory is provided by the closed- 
time-path (CTP) method. Here all quantities are given on a special real-time contour 



C(x) = -d^(x)d^(x) 




A 
4! 




(1) 
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with the time argument running from — oo to +00 on the chronological branch (+) and 
returning from +00 to —00 on the anti-chronological branch (— ). In particular the path- 
ordered Green functions are defined as {x\ = (ti,xi), x 2 = (t 2 ,x 2 )) 



G(xi,x 2 ) = (TP{ 4>{x x ) <f>(x 2 )}) (2) 
= Q P (ti - h) ( <j)(x 1 ) <p(x 2 ) ) + e p (t 2 - h) ( cj)(x 2 ) <p(xi) ) , 

where the operator T p orders the field operators according to the position of their argu- 
ments on the real-time path as accomplished by the path step-functions P 

1 , if ti is later on the contour as t 2 , 

e p (ti-t 2 ) = { (3) 

, if t 2 is earlier on the contour as t% . 

The expectation value in (J2J) is taken with respect to some given density matrix po, which 
is constant in time, while the operators in the Heisenberg picture contain the whole 
information of the time dependence of the nonequilibrium system. 

Self-consistent equations of motion for these Green functions can be obtained with 
help of the two-particle- irreducible 1 (2PI) effective action T[G], which is given by the 
double Legendre transform of the generating functional of connected Green functions W 
with respect to the local and bilinear external sources |51j . The functional which is 
the sum of all closed 2PI diagrams built up by full propagators G (and bare vertices) 
|2~%t ESI E2j , determines the resummed self-energy by functional variation as 

E( ^ 2) - % m^D- (4) 

In general, the self-energy E contains a singular contribution S 5 on the contour and a 
non-local part, which can be expressed by a sum over 9 p -functions 

E(x!,x 2 ) = J: 5 (x 1 )5lf +1 \x 1 -x 2 ) (5) 
+ Q p {ti-t 2 ) E > (x 1 ,x 2 ) + e p (t 2 -*i) S < (xi,x 2 ) . 

Here the <5 p -function on the closed-time-path is defined as 

5f +l \ Xl - x 2 ) = 6 p (h - t 2 ) 5^{x x - x 2 ) , (6) 

{$(t\ — h) , if both t\ and t 2 are on the upper (+) branch , 
—5(ti — t 2 ) , if both t\ and t 2 are on the lower (— ) branch , 
, in other cases , 



^-particle-irreducible diagrams remain connected when n arbitrary internal lines are cut. 
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in accordance with the contour integration while d denotes the spatial dimension of the 
problem. The functional $ allows - by restricting to a particular set of diagrams - for the 
construction of effective theories, which are conserving and thermodynamically consistent 
15%] . In our present calculations we will take into account contributions up to the 
three-loop order for the ^-functional (cf. Fig. 1 in Ref. |51j). Whereas the self-energies 
are defined by (jlj) the evolution of the Green function is determined by the stationarity 
condition 5Y/5G = jHHEO- 

For actual calculations it is advantageous to change to a single-time representation for 
the Green functions and (the non-local part of) the self-energies defined on the closed- 
time-path. In line with the position of the coordinates on the contour there exist four 
different two-point functions 



i G c (xi, x 2 ) 
i G < (xi, £2) 
i G > (xi, X2) 
iG a (x 1 ,x 2 ) 



i G ++ (xi, X2) 
i G + ~(xi, X2) 
i G (xi, X2) 
i G (xi, x 2 ) 



1 (sr i^X ^ j 37 2 

1 Cj i^X -y j 37 2 

1 Cj 1 J 3^2 

X Cj 1 J 3^2 



(T c {<t>{x x )j>{x 2 )}), 
{{(f){x 2 ) 0(Xi)}), 

({0(Xi)0(x 2 )}), 
(^{0(^)^2)}). 



(7) 



Here T c (T a ) represent the ( ant i-) time-ordering operators in case of both arguments lying 
on the (anti-) chronological branch of the real-time contour. These four functions are not 
independent of each other. In particular the non-continuous functions G c and G a are 
built up by the Wightman functions G > and G < and the usual O-functions in the time 
coordinates as 



G%x 1 ,x 2 ) = Q{t 1 -t 2 )G>(x 1 ,x 2 ) + Q{t 2 -t 1 )G<(x 1 ,x 2 ), (8) 
G a {x u x 2 ) = e(t 2 -ti) G > (x 1 ,x 2 ) + 6(ti-t 2 ) G<(x u x 2 ). 



By fixing the time arguments of the Green functions on different branches of the 
contour and resolving the time structure of the path-ordered quantities we obtain the 
Kadanoff-Baym equations for the time evolution of the Wightman functions I51| 
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- [dpd^+m 2 ] G%( Xl ,x 2 ) = ^ 5 (x 1 )G^(x 1 ,x 2 ) 



(9) 



+ dz d d z [E > (xi,z)-E < (xi,z)\ G^(z,x 2 ) 



- J dzo Jd d z t£(x!,z) [G > (z,x 2 )-G < (z,x 2 ] 



- [d£ 2 d% 2 + m 2 ] G*( Xl ,x 2 ) = Z?(x 2 )G*(xi,x 2 ) 

+ JdzQ jd d z [G > {x 1 ,z)-G < (x 1 ,z)} E^(z,x 2 ) 



'to 

fdz !d d z G*{xi,z) [E > {z,x 2 ) - E<(z,x 2 ) 

J t J 



Within the three-loop approximation for the 2PI effective action $ we get two different 
self-energies: In leading order of the coupling constant only the local tadpole diagram 
contributes and leads to the generation of an effective mass for the field quanta. This 
self-energy (in coordinate space) is given by 

= ^tG < (x,x). (10) 

In next order in the coupling constant (i.e. A 2 ) the non-local sunset self-energy (cf. Fig. 
2 in jHT]) enters the time evolution as 

Y£{x u x 2 ) = -^G*(x 1 ,x 2 )GZ(x 1 ,x 2 )G$(x 2 ,x 1 ) = ~ [G*(x lt x 2 )] a . (11) 





Thus the Kadanoff-Baym equation (jHJ) in our case includes the influence of a mean-field 
on the particle propagation - generated by the tadpole diagram - as well as scattering 
processes as inherent in the sunset diagram. 

The Kadanoff-Baym equation @ describes the full quantum nonequilibrium time 
evolution on the two-point level for a system prepared at an initial time to, i.e. when 
higher order correlations are discarded. The causal structure of this initial value problem 
is obvious since the time integrations are performed over the past up to the actual time 
t\ (or t 2 , respectively) and do not extend to the future. 

Furthermore, also linear combinations of the Green functions in single-time represen- 
tation are of interest and will be exploited for the spectral properties of the system later 
on. The retarded Green function G R and the advanced Green function G A are given as 
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G R (x 1 ,x 2 ) = Q(h-t 2 ) [G > (x 1 ,x 2 )-G < (x 1 ,x 2 )] (12) 

= e(ti-t 2 ) ([0(xi), 4>{x 2 )\_) 

= G c (xi,x 2 ) - G < (x 1 ,x 2 ) = G > (x 1 ,x 2 ) - G a (x 1 ,x 2 ), 

G A ( Xl ,x 2 ) = -Qfo-h) [G>(x 1 ,x 2 )-G<(x 1 ,x 2 )} (13) 
= -e(t 2 -ti) ([tixx) , <f>(x 2 ))_) 

= G c (x\,x 2 ) — G > (xi,x 2 ) = G < (xi,x 2 ) — G a (xi,x 2 ). 

These Green functions contain exclusively spectral, but no statistical information of the 
system. Their time evolution is determined by Dyson- Schwinger equations and given by 
(cf. Ref. jSU) 

- [^ 1 ^ l + m 2 + S <5 (x 1 )j G R ( Xl ,x 2 ) = (14) 

5( d+1 \ Xl -x 2 ) + [d d+1 z E R (x h z) G R (z,x 2 ), 



~ [d^d^ 1 +m 2 + E s {x 1 )] G A ( Xl ,x 2 ) = (15) 

5 (d+1) (x!-x 2 ) + jd d+l z Z A (x u z) G A (z,x 2 ), 

where the retarded and advanced self-energies J2 R , T, A are defined via S > , S < similar to 
the Green functions ()12j) and ([13)1 . Thus the retarded (advanced) Green functions are 
determined by retarded (advanced) quantities, only 

2.2 Generalized transport equations 

The derivation of generalized transport equations has been presented before in Refs. [21 
El EEH 1221 ESI ESI Ell EHj- We here briefly repeat the derivation in order to introduce 
the necessary notations and to work out the differences in the Kadanoff-Baym 9\ and 
Botermans-Malfliet [22] scheme that will be investigated in Section 3. 

We start by rewriting the Kadanoff-Baym equation for the Wightman functions in 
coordinate space (xi = (t 1 ,xi),x 2 = (t 2 ,x 2 )) (jHJ) as 

[d^d* 1 + m 2 + T 1 s {x l )\ iG^( Xl ,x 2 ) = ilf( Xl ,x 2 ). (16) 

The collision terms on the r.h.s. of (|16j) are given in D = d + 1 space-time dimensions by 
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convolution integrals over coordinate space self-energies and Green functions: 



I?(x 1} x 2 ) = f^z [^>(x l ,z)-^ < (x l ,z)\ G^(z,x 2 ) (17) 
+ fd D z T^{x u z) [G > {z,x 2 )-G < {z,x 2 )} . 



In the general case of an arbitrary (scalar) quantum field theory S 5 is the local (non- 
dissipative) part of the path self-energy while £< resemble the non-local collisional self- 
energy contributions. In the representation (jl7j) the integration boundaries are exclusively 
given for the time coordinates since the integration over the spatial coordinates extends 
over the whole spatial volume from — oo to +00. 

Since transport theories are formulated in phase-space one changes to the Wigner 
representation via Fourier transformation with respect to the rapidly varying ('intrinsic') 
relative coordinate Ax = x\ — x 2 and treats the system evolution in terms of the ('macro- 
scopic') mean space-time coordinate x = (x\ +x 2 )/2 and the four-momentum p = (po,p). 
The functions in Wigner space are obtained as 

/oo 
d D Ax e + iA W F(x! =x + Ax/2, x 2 = x- Ax/2). (18) 
-00 

For the formulation of transport theory in the Wigner representation we have to focus not 
only on the transformation properties of ordinary two-point functions as given in (|18j) . 
but also of convolution integrals as appearing in Eq. (|17|). A convolution integral in D 
dimensions (for arbitrary functions F, G) , 

/oo 
d D z F(x h z) G(z,x 2 ) (19) 
-00 



transforms as 



H(p, x) 



/oo 

d D Ax e +iAx » p " H(x 1 ,x 2 ) (20) 
-OO 

/oo roc 

d D Ax e +iAx » p " / d D z F{x u z) G(z,x 2 ) 
-oo J —00 



[ F(p,x) G(p',x ) ] 



In accordance with the standard assumption of transport theory we assume that all func- 
tions only smoothly evolve in the mean space-time coordinates and thus restrict to first 
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order derivatives. All terms proportional to second or higher order derivatives in the mean 
space-time coordinates (also mixed ones) will be dropped. Thus the Wigner transformed 
convolution integrals (|19|) are given in first order gradient approximation by, 

H(p,x) = F(p,x) G(p,x) + i^{F(p,x),G(p,x)} + 0{d 2 x ) , (21) 
using the relativistic generalization of the Poisson bracket 

{ F(p, x) , G(p, x) } = F(p, x) ■ <9£ G(p, x) - % F(p, x) ■ ^ G(p, x) . (22) 



In order to obtain the dynamics for the spectral functions within the approximate 
scheme we start with the Dyson-Schwinger equations for the retarded and advanced Green 
functions in coordinate space (|14I15|) . - We note that the convolution integrals in (|14|) and 
(|15|) extend over the whole space and time range in contrast to the equations of motion 
for the Wightman functions given in (jTS)) and (|17j). - The further procedure consists in 
the following steps: First we 

i) transform the above equations into the Wigner representation and apply the first order 
gradient approximation. In this limit the convolution integrals yield the product terms and 
the general Poisson bracket of the self-energies and the Green functions { H R ' A , G R ^ A }. 
We, furtheron, represent both equations in terms of real quantities by the decomposition 
of the retarded and advanced Green functions and self-energies as 

qR/a = Re QR ± lIm QR = Re QR T iAjl, A = T2ImG R / A , 

(23) 

e r/a _ Re Y, R ± i/mE R = ReY* R iV/2 , T = ^2JmS R / A . 

We find that in Wigner space the real parts of the retarded and advanced Green functions 
and self-energies are equal, while the imaginary parts have opposite sign and are propor- 
tional to the spectral function A and the width f, respectively. The next step consists in 

ii) the separation of the real part and the imaginary part of the two equations for the 
retarded and advanced Green functions, that have to be fulfilled independently. Thus 
we obtain four real- valued equations for the self-consistent retarded and advanced Green 
functions. In the last step 

iii) we get simple relations by linear combination of these equations, i.e. by adding/subtrac- 
ting the relevant equations. 

This finally leads to two algebraic relations for the spectral function A and the real 
part of the retarded Green function ReG R in terms of the width F and the real part of 
the retarded self-energy i?eS fl as 
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[p 2 -p 2 -m 2 -Z s + ReZ R ] ReG R = 1 + ~T A, (24) 
[Po-P 2 ~m 2 -t, s + ReZ R ]A = f ReG R . (25) 

Note that all terms with first order gradients have disappeared in ()24|) and (|25|). A first 
consequence of f!25|) is a direct relation between the real and the imaginary parts of the 
retarded/advanced Green function, which reads (for f ^ 0): 

KG* A. (26) 

Inserting (}26|) in (|2*4*j) we end up with the following result for the spectral function and 
the real part of the retarded Green function 



A [p 2 - p 2 - m 2 - Z s -i?eS«] 2 + r 2 /4 M 2 + r 2 /4' ^ 

R [p 2 ~p 2 -m 2 -f: s -ReJ: R ] M 

[p 2 -p 2 -m 2 -Z s -i?eS R ] 2 + f 2 /4 ~ M 2 + f 2 /4' 1 ' 

where we have introduced the mass-function M in Wigner space: 



M(p,x) = p 2 -p 2 -m 2 -X 5 (x)-ReZ R (p,x) . (29) 

The spectral function (}2Tj) shows a typical Breit- Wigner shape with energy- and momentum- 
dependent self-energy terms. Although the above equations are purely algebraic solutions 
and contain no derivative terms, they are valid up to the first order in the gradients. 

In addition, subtraction of the real parts and adding up the imaginary parts lead to 
the time evolution equations 

p»d%A = l -{Y> 8 + Re^ R , A} + i{f, ReG R } } (30) 

p^d^ReG R = ]-{Z s + ReZ R , ReG R } - i{f , A}. (31) 

The Poisson bracket containing the mass-function M leads to the well-known drift oper- 
ator p^ d x F (for an arbitrary function F), i.e. 
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{M , F} = {pl-f -rn 2 -E 5 - ReZ R , F} (32) 
= 2jfd^F — {E s + ReE R , F} , (33) 



such that the first order equations (fHUJ) and (|3*T|) can be written in a more comprehensive 
form as 



{M , A} = {T , Re G R } , (34) 

{ M , ReG R } = -I{r,A}. (35) 

When inserting ()27)1 and ()28j) we find that these first order time evolution equations are 
solved by the algebraic expressions. In this case the following relations hold: 

M 2 -f 2 /4 

{M,A} = {F,ReG R } = {M,r} [WTf ^ , (36) 

{M , Re G } = -j{r,i} = {M,r} [A?2 + j 4 /4|2 . (37) 

Thus we have derived the proper structure of the spectral function (j27J) within the first- 
order gradient (or semiclassical) approximation. Together with the explicit form for the 
real part of the retarded Green function ()28)l we now have fixed the dynamics of the 
spectral properties, which is consistent up the first order in the gradients. 

As a next step we rewrite the memory terms in the collision integrals ()17|) such that 
the time integrations extend from — oo to +oo. In this respect we consider the initial time 
to = — oo whereas the upper time boundaries ti,t 2 are taken into account by G-functions, 
i.e. 

/oo 
d D x' 6(ti -t') [E > (xi,x') - S < (x 1 ,x / )] G^(x',x 2 ) 
-oo 



■oo 



+ / d D x' £^(xi,x') e(t 2 -t') [G > (x',x 2 )-G < (x',x 2 



/oo 
d D x' E R {x u xf) G^{x',x 2 ) + Y^{ Xl ,x') G A (x',x 2 ) . (38) 
-oo 
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We now perform the analogous steps as invoked before for the retarded and advanced 
Dyson- Schwinger equations. We start with a first order gradient expansion of the Wigner 
transformed Kadanoff-Baym equation using (|38|) for the memory integrals. Again we sep- 
arate the real and the imaginary parts in the resulting equation, which have to be satisfied 
independently. At the end of this procedure we obtain a generalized transport equation 

2p»d£ iG^ - {Z 5 +ReZ R , iG^} - , ReG R } = zS< iG> - iT? iG< 



{M,iG%} - {iY^, ReG R } = iY< iG> - zE> %G< (39) 

as well as a generalized mass-shell equation 

[p 2 -m 2 -E s - Re£ R ] iG* = if^ Re G R + J{*£> iG< } - ^{iS<, iG> } (40) 



with the mass-function M specified in (|2*5|) . Since the Green function G^(xi,x 2 ) consists 
of an antisymmetric real part and a symmetric imaginary part with respect to the relative 
coordinate X\ — x 2 , the Wigner transform of this function is purely imaginary. It is thus 
convenient to represent the Wightman functions in Wigner space furtheron by the real- 
valued quantities iG<(p,x). Since the collisional self-energies obey the same symmetry 
relations in coordinate space and in phase-space, they will be kept also as ill<(p, x) 
furtheron. 

In the transport equation (}3*9"j) one recognizes on the l.h.s. the drift term p^ iG<, 
as well as the Vlasov term with the local self-energy and the real part of the retarded 
self-energy Re Y, R . On the other hand the r.h.s. represents the collision term with its 
typical 'gain and loss' structure. The loss term iS > iG < (proportional to the Green 
function itself) describes the scattering out of a respective phase-space cell whereas the 
gain term iS < iG > takes into account scatterings into the actual cell. The last term on 
the l.h.s. { iT^, Re G R } is very peculiar since it does not contain directly the distribution 
function %G < . This second Poisson bracket vanishes in the quasiparticle approximation 
and thus does not appear in the on-shell Boltzmann limit. As demonstrated in detail in 
Refs. El M EH EZ1 EE1 the second Poisson bracket { zE^, Re G R } governs the evolution 
of the off-shell dynamics for nonequilibrium systems. 

Although the generalized transport equation (J39|) and the generalized mass-shell equa- 
tion (|4U|) have been derived from the same Kadanoff-Baym equation in a first order gradi- 
ent expansion, both equations are not exactly equivalent 123 EH1 • Instead, they deviate 
from each other by contributions of higher gradient order (see below and Refs. 122j). 
This raises the question: which one of these two equations has be considered to higher 
priority? The question is answered in practical applications by the prescription of solv- 
ing the generalized transport equation (|39|) for iG < in order to study the dynamics of 
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the nonequilibrium system in phase-space. Since the dynamical evolution of the spectral 
properties is taken into account by the equations derived in first order gradient expansion 
from the retarded and advanced Dyson-Schwinger equations, one can neglect the general- 
ized mass-shell equation (J4Up . Thus for our actual numerical studies in Section 3 we will 
use the generalized transport equation (|3T?|) supported by the algebraic relations (}2*7|) and 



Furthermore, one recognizes by subtraction of the iG > and iG < mass-shell and trans- 
port equations, that the dynamics of the spectral function A = iG > — iG < is determined 
in the same way as derived from the retarded and advanced Dyson-Schwinger equations 
(J2ZJ) and (jHIJ). The inconsistency between the two equations (fBTiJ) and (|4~Uj) vanishes since 
the differences are contained in the collisional contributions on the r.h.s. of (|39p. 

In order to evaluate the { iS < , Re G R }-term on the l.h.s. of (}3*9"j) it is useful to intro- 
duce distribution functions for the Green functions and self-energies as 



iG < (p,x) = N(p,x) A(p,x) 



iG > (p,x) = [1 + N(p,x)] A(p,x) 



(41) 



iE < (p, x) = iV E (p, x) x) 



iS>(p,x) = [1 + N s (p,x)} T(p,x) . (42) 



In equilibrium the distribution function with respect to the Green functions N and the 
self-energies iV E are given as Bose functions in the energy p at given temperature and 
thus are equal. Following the argumentation of Botermans and Malfliet |22] the distribu- 
tion functions N and iV s in (141 J) should be identical within the second term of the l.h.s. 
of ()39|) within a consistent first order gradient expansion. In order to demonstrate their 
argument we write 



ill' 



r n 1 



r n + k 



(43) 



The 'correction' term 



K 



r ( n - N' 



iE< iG> - iE> iG< ) A' 



(44) 



is proportional to the collision term of the generalized transport equation (}3*9"j) . which 
itself is of first order in the gradients. Thus, whenever a distribution function iV s appears 
within a Poisson bracket, the difference term (iV 2 — N) becomes of second order in the 
gradients and should be omitted for consistency. As a consequence iV s can be replaced 
by N and thus the self-energy E < by G < ■ T/A in the Poisson bracket term {E < , ReG R }. 
The generalized transport equation (}3*9"j) then can be written in short-hand notation 



l -Af 
2 



{ M, iG<} - i { T , M • iG < } 



iE< iG> - zE> iG' 



(45) 
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with the mass-function M (|29p . The transport equation PB^) within the Botermans- 
Malfliet (BM) form resolves the discrepancy between the generalized mass-shell equation 
(|4U|) and the generalized transport equation in its original Kadanoff-Baym (KB) form 
flUJ). 

We recall that the consistent Botermans-Malfliet form ()45|) has been used as a starting 
point for the off-shell description of realistic heavy-ion collisions (21 El HH] • Together with 
an extended test-particle ansatz for iG < this equation allows for a convenient solution of 
the generalized transport equation in phase-space including the dynamics of the spectral 
function. However, it is presently not clear in how far the generalized transport equa- 
tions in KB-form or in BM-form (|43|l reproduce the same dynamics as for the full 
Kadanoff-Baym theory Moreover, the differences in the time evolution of nonequi- 
librium systems between the KB-form (JHHJ) and BM-form ()45j) are not known either. We 
will thus perform an explicit comparison between the different limits for 4 -theory in 2+1 
dimensions in the following Section. For a related and more detailed comparison between 
the solution of the full KB equations with those from the on-shell Boltzmann limit we 
refer the reader to Section V of Ref . [5T] . 

3 Dynamics within the Generalized Transport Equa- 
tion 

In this Section we will perform numerical studies of the dynamics inherent in the gen- 
eralized transport equations derived in Section 2.2 in comparison to solutions of the full 
Kadanoff-Baym equations (jHJ) for the </> 4 -theory in 2+1 space-time dimensions within the 
three-loop approximation for the effective action carried out in Ref. [ST]. This fixes the 
self-energies E" 5 and iT^ in flU)) and (gUJ) to be the same as in the case of the full Kadanoff- 
Baym theory given by the tadpole (JTUJ) and the sunset (fTTj) contributions, respectively. 
We recall that both self-energies are ultraviolet divergent and have to be renormalized by 
proper counter terms as shown in Appendix B of Ref. [5TJ. In the following we will use 
the renormalized mass m—1, which implies that times are given in units of the inverse 
mass. Accordingly, the coupling A is given in units of mass m such that tm and X/m are 
dimensionless. 

For the first investigation we concentrate on the dynamics of the generalized trans- 
port equation As in Ref. [STJ we restrict ourselves to homogeneous systems in space. 
Consequently the derivatives with respect to the mean spatial coordinate x vanish, such 
that the generalized transport equation reduces to 



2p d t iG < - {E s +ReE R ,iG K } - {i£< , ReG R } = iS< iG > - iT> iG< (46) 



with the simplified Poisson brackets (for arbitrary functions F, G) 
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{F(p,po,t) , G(p,po,t) } 



(47) 



= dp Q F(P, Po, t) d t G(p, p , t) - d t F(p, p , t) d VQ G(p, p Q ,t) . 

In order to obtain a numerical solution the KadanofT-Baym equation © as well as the 
generalized transport equations are transformed to momentum space. For the numer- 
ical solution we have developed a flexible and accurate algorithm, which is described in 
more detail in Appendix A of Ref. |51| . Furthermore, since in 2+1 space-time dimensions 
both self-energies (jlU|) . (jll|) are ultraviolet divergent, they have to be renormalized by in- 
troducing proper counter terms. The details of the renormalization scheme are given in 
Appendix B of Ref. [5T] as well as a numerical proof for the convergence in the ultraviolet 
regime. For the actual numerical integration of the generalized transport equation we use 
the same grid in momentum space as employed for the full KB-theory in Ref. [ST] • 

3.1 The Initial State 

In order to specify the problem for the first order transport equation in time, the initial 
state has to be fixed, i.e. the initial Green function iG < (p,p , t = 0) has to be specified for 
all momenta p and all energies po- To this aim, we first assume a given initial distribution 
function n(p, t = 0) in momentum space equivalent to those used in the investigation of 
the full Kadanoff-Baym theory in Ref. [5T] in order to allow for a comparison of the two 
schemes. Next, the complete initial phase-space distribution function N(p,p ,t = 0) has 
to be specified as a function of the energy p as well. We achieve this by assuming N 
to be constant in the energy variable po as a starting point. In order to determine the 
initial Green function we employ the self-consistent iteration procedure for the full spec- 
tral function as described in Appendix D of Ref. jST] using the specified nonequilibrium 
distribution function N. The iteration process then yields the fully self-consistent spectral 
function A(p, p ,t = 0) for this initial distribution and thus determines the initial Green 
function via 

iG < (p,p ,t = 0) = N(p,p ,t=0) A(p,p ,t=0). (48) 

In principle one might choose the initialization - in particular the energy dependence 
of the distribution functions - in various ways. We adopt here the simple description 
introduced above since it approximately reproduces the same equal-time (i.e. energy 
integrated) initial Green function as in the full Kadanoff-Baym calculation in Ref. [oT] . 
However, with increasing coupling strength A the spectral function achieves a considerable 
width and the on-shell energies of the momentum modes are slightly shifted. In this case 
the equal-time Green function deviates somewhat in the low momentum region from the 
'quasiparticle'-like initialization used for the Kadanoff-Baym calculation (cf. Ref. [oT]). 

The initial distribution functions for the following studies are shown in Fig. Q for 
polar symmetric configurations as a function of the momentum coordinate p x for p y = 0. 
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They are identical to those used in Ref. jSTj within the full Kadanoff-Baym theory. The 
resulting initial equal-time Green functions - for the coupling constant X/m — 18 - are 
displayed in Fig. |21 By comparison with the initial Green functions of the Kadanoff-Baym 
calculation (c.f. lower part of Fig. 3 in Ref. [HI]) we find tiny deviations in the region of 
small absolute momenta. This is, of course, a consequence of the large coupling constant 
A employed, which we use in order to compare to the calculations within the full quantum 
evolution for the same coupling strength. Since the initial states are very close to those 
used for the Kadanoff-Baym theory, we will also denote them as initializations Dl, D2 
and D 3 as in Ref. [5T] . 

The advantage of the initialization prescription introduced above is that the actual 
spectral function - directly obtained by A = iG > — iG < from the Green functions - com- 
plies with the one determined from the self-energies (|27|) in accordance with the first order 
gradient expansion scheme. During the nonequilibrium time evolution this correspondence 
is maintained since the analytic expression for the spectral function already is a solution 
of the generalized transport equation itself. Furthermore, the real part of the retarded 
Green function, that enters the peculiar second Poisson bracket on the l.h.s. in ()39j) . can 
be taken in the first order scheme ()28|) which simplifies the calculations considerably. 

We mention that other prescriptions are also possible for the calculation of Re G R : 

1) One can determine it directly from iG^ by a Fourier technique similar to one used in 
the calculation of the self-consistent spectral functions in Appendix D of Ref. |5T]- Here 
the real part of the retarded Green function is obtained via inverse Wigner transformation 
with respect to the energy, multiplication by the B-function in relative time and trans- 
formation back to phase-space. 

2) One can also use dispersion relations with the spectral function to specify Re G R . 

As we have checked in our actual simulations all prescriptions lead to practically identical 
results. 

3.2 Numerical Study of Equilibration 

Now we turn to the actual solutions of the generalized transport equation in the KB form 
(|39|). In Fig. El (upper part) we show the time evolution of the equal-time Green function 
iG < (\p\,t,t) for the polar symmetric initial states Dl, D2 and D3 as specified in Fig. 
2. We have displayed several momentum modes \p\/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 of 
the equal-time Green function on a logarithmic time scale. As in the full Kadanoff-Baym 
theory (lower part) we find that for all initializations the quantum system approaches a 
stationary state for t — > oo, i.e. all momentum modes approach a constant. 

However, the respective momentum modes of the different initializations do not achieve 
identical values for t — » oo, as seen in particular for the low momenta \ p\/m = 0.0, 0.8 
in Fig. El (upper part). This is not surprising since the various initializations - obtained 
within the self-consistent scheme described above - do not correspond to exactly the same 
energy. This is why the respective long-time limits differ slightly. The small difference 
in energy is, of course, most prominently seen in the low momentum (energy) modes. 
Moreover, the dynamics within the generalized transport equation ()39|) is in general very 
similar to the full Kadanoff-Baym theory (lower part). For all three initial states we find 
(apart from the very initial phase t ■ m < 5) the same structures during the equilibration 
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process. In particular for the initializations Dl and D2 the characteristic overshooting for 
the low momentum modes is seen as in the full quantum evolution, which does not show 
up in solutions of the corresponding on-shell Boltzmann limit. Since in the Boltzmann 
limit a strictly monotonous evolution of the momentum modes is seen (cf. Ref. j51j ) 
this overshooting has to be attributed to an off-shell quantum effect. Even the positions 
of the maxima are in a comparable range: For the initialization Dl they are shifted to 
slightly larger times and are a little bit lower than in the full calculation; the same holds 
for the initial state D2. The initial distribution D3 yields a monotonous behaviour for all 
momentum modes within the generalized transport formulation which is again in a good 
agreement with the full dynamics. 

Some comments are worthwhile with respect to the comparison performed above: The 
spectral function in the Kadanoff-Baym calculation is completely undetermined in the 
initial state; it develops during the very early phase to an approximate form (which in the 
following still evolves in time). In contrast to this, the spectral function in the generalized 
transport formulation (|39|) has a well-defined structure already from the beginning. This 
principle difference results from the fact, that in the Kadanoff-Baym case we deal with 
a true initial value problem in the two time directions (ti,t 2 ). Moreover, the spectral 
distribution for very early times depends on the specific initial conditions adopted (cf. 
Appendix C of Ref. |51j). Additionally, the relative time integral in (ti— £2) ~~ to obtain the 
spectral function in energy po by Wigner transformation - is very small in the initial phase 
as discussed in Ref. [HI] • Consequently, the spectral shape in Wigner space is determined 
by the finite integration interval in time rather than by the interactions itself. On the other 
hand, we have used an infinite relative time range in deriving the generalized transport 
equation within the first order gradient expansion. Thus in this case we deal with a 
completely resolved spectral function already at the initial time, that clearly exhibits the 
physics incorporated. This demonstrates why both approaches can only be compared to 
a certain extent for the very early times. 

Nevertheless, the inclusion of the dynamical spectral function in terms of the general- 
ized transport equation surpasses the shortcomings of the quasiparticle Boltzmann limit 
discussed in Section V and Appendix E of Ref. |5T]- Whereas the latter approach leads 
to a strictly monotonous evolution of the momentum modes, the inclusion of quantum 
effects in terms of a semiclassical approximation correctly yields the overshooting effects 
as in the solution of the full Kadanoff-Baym equation. 

Finally, concentrating on the very early time behaviour, we find a significant difference 
between the full and the approximate dynamics in the gradient scheme (|39|) . For the 
generalized transport equation we find a monotonous evolution of the equal-time Green 
function momentum modes, whereas strong oscillations are observed in the initial phase 
for the solution of the full Kadanoff-Baym theory (see also Fig. 4 of Ref. |51j). Thus, with 
respect to the early time behaviour the generalized transport equation behaves much more 
like the Boltzmann approximation, which is a first order differential equation in time as 
well. However, the Kadanoff-Baym evolution is given by an integro-differential equation 
of second order in time. In this case the phase correlations between the Green functions 
G^, GS, G< n are kept and the instantaneous switching-on of the interaction results in 
an oscillatory behaviour of the single momentum modes. We mention that correlations 
build up very rapidly in the very early phase of the evolution in the full KB theory 
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(cf. Fig. 6 of Ref. [3T]) that may be attributed to this oscillatory behaviour. When 
including the collisional self-energies, i.e. on the three-loop level for the effective action, 
these oscillations are damped in time typically with the respective on-shell width for given 
momentum. If only the tadpole term is included, these oscillations are not damped and 
maintain forever (see also Fig. 4 of Ref. [51 J ) . Therefore, the origin of the oscillatory 
behaviour can be traced back to the order of the underlying differential equation and the 
rapid build-up of correlations that have not been incorporated in the initial conditions for 
the full KB-theory. 

3.3 Evolution of the Spectral Function 

Since the Green functions develop in time also the spectral properties of the system change 
as well. In Fig. EJthe time evolution of the spectral functions for the initializations Dl, 
D2 and D3 within the gradient scheme are displayed. We focus on the spectral functions 
for two particular momentum modes \p\/m = 0.0 (l.h.s.) and \p\/m = 2.0 (r.h.s.) for 
various system times t ■ m — 5, 20, 60, 120, 240, 360, 480, 600 up to the long-time limit. 
This representation corresponds to Fig. 7 in Ref. |51j . where the respective evolution of 
the spectral function is studied for the full Kadanoff-Baym theory. We find that the time 
evolution of the spectral functions obtained from the generalized transport equation (|39|) 
is very similar to the one from the full quantum calculation (see below). The zero- mode 
spectral function for the initial distribution Dl becomes sharper with time and is moving 
to slightly higher energies. The opposite characteristics is observed for the zero-mode 
spectral function for the initialization D3, which broadens with time (reducing the peak 
correspondingly) and slowly shifts to smaller energies. Together with the weak evolution 
for the distribution D2 (which only slightly broadens at intermediate times and returns 
to a narrower shape at smaller energies in the long-time limit) the evolution of all three 
initializations in the semiclassical approximation is well comparable to the full Kadanoff- 
Baym dynamics (cf. Fig. 7 in Ref. |3T]). Furthermore, the maxima of the zero-mode 
spectral functions are located above the bare mass (as indicated by the on-shell arrow) 
for all initial states during the time evolution. 

The spectral functions for the momentum mode \p\jm = 2.0 are in a good agreement 
with the Kadanoff-Baym dynamics as well. Again we observe - for the initial distribution 
Dl - a narrowing of the spectral function, while for D3 the spectral function broadens 
with time. Moreover, the width of the spectral function starting from distribution D2 
shows a non-monotonous behaviour with a maximum at intermediate times. 

In order to study the dynamics of the spectral function in a more quantitative manner 
we display in Fig. El the time evolution of the on-shell energies (as derived from the 
maxima of the spectral function) for the momentum modes \p\/m = 0.0 (upper plot) 
and 2.0 (lower plot) for the initializations Dl, D2 and D3 with X/m = 18 (l.h.s.). By 
comparison with the corresponding results from the Kadanoff-Baym theory (r.h.s.) we 
observe a close similarity of the evolutions within the full and the semiclassical KB scheme. 
The effective mass of the zero momentum mode decreases for initialization D3, passes a 
minimum for D2 and increases for the initial state Dl. 

As familiar from the Kadanoff-Baym calculations in Section IV. B of (51] the behaviour 
of the on-shell energies is different for higher momentum modes. We find for the momen- 
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turn mode \p\jm = 2.0 a monotonous decrease of the on-shell energy for the initializations 
Dl and D2 and an increase for distribution D3. Altogether, the evolution of the on-shell 
energies for the higher modes is rather moderate compared to the lower ones in accordance 
with the dominant momentum contribution and the weakening of the retarded self-energy 
for higher energy modes. 

Finally, the on-shell energies approach a stationary state for all modes and all initial- 
izations. However, the long-time limit of the equal momentum modes is not exactly the 
same for all initial distributions Dl, D2 and D3. As discussed above this small differ- 
ence can be traced back to the specific initial state generation from the given momentum 
distribution. 

Next we consider the time evolution of the on-shell width as determined by the imag- 
inary part of the retarded self-energy at the maximum position of the spectral function. 
In Fig. |U1 (l.h.s.) the on-shell width is displayed for the two momentum modes | p\/m = 
0.0 and \p\/m = 2.0 for all three initial distributions Dl, D2 and D3 with X/m =18 
as a function of time. For both momentum modes the on-shell width increases for the 
distribution D3, while it has a maximum at intermediate times (t ■ m « 40) for the ini- 
tialization D2. Thus the results - together with the reduction of the on-shell width for 
both momentum modes for the initialization Dl - is in good agreement with the results 
obtained for the full Kadanoff-Baym theory (r.h.s.). However, the stationary values for 
the on-shell widths deviate again slightly in accordance with the different preparation of 
the initial state in the gradient scheme. 

In summarizing we find that the main characteristics of the full quantum evolution of 
the spectral function are maintained in the semiclassical transport equation ()39|) as well. 
This includes the evolution of the on-shell energies as well as the width of the spectral 
function. Since the generalized transport equation is formulated directly in Wigner space 
one has access to the spectral properties at all times, whereas the very early times in the 
Kadanoff-Baym case have to be excluded due to the very limited support in the relative 
time interval (t\ — £2) for the Wigner transformation. 

3.4 Stationary State of the Semiclassical Evolution 

As we have observed in the previous Subsections the evolution within the generalized 
transport equation (}39|) leads to a stationary state for all three different initializations 
Dl, D2 and D3. Thus we turn to the investigation of this long-time limit itself, here in 
particular for the initialization D2. In Fig. (l.h.s.) we show the distribution function N 
of various momentum modes \p\/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 for large times (t ■ m = 
600) as derived from the Green function itself and the spectral function via the relation 
N = iG < I A. The distribution function for a given momentum mode is calculated for 
all energies Pq where the corresponding spectral function - as displayed in the lower 
part of Fig. [7| exceeds a value of 0.5. Since the width of the late time spectral function 
decreases with increasing momentum, the energy range for which the distribution function 
is shown, is smaller for larger momentum modes. We find, that all momentum modes of 
N can be fitted at all energies by a single Bose function with a temperature T/m = 1.827 
and a very small chemical potential fi/m = 0.000817. Thus the generalized transport 
formulation (|39|) leads to a complete (off-shell) equilibration of the system very similar to 
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the solution of the full Kadanoff-Baym equation (r.h.s. of Fig. |7J). Furthermore, the long- 
time limit of the semiclassical time evolution exhibits a vanishing chemical potential fi/m 
in accordance with the properties of the neutral 4 -theory. This might have been expected 
since in the generalized transport equation particle number non- conserving processes of 
the type 1 <-> 3 - which lead to the decrease of the chemical potential - are included 
by means of the dynamical spectral function. Thus the semiclassical approximation (}3T?j) 
solves the problems within the Boltzmann limit, which does not yield a relaxation of the 
chemical potential, since only on-shell 2^2 transitions of quasiparticles are taken into 
account as demonstrated in Section V of Ref . [HI] • 

After observing, that the chemical potential decreases to zero in the long-time limit, 
it is interesting to study the relaxation process itself. The relaxation of the chemical 
potential fi/m is shown for the three different initializations Dl, D2 and D3 with coupling 
constant A/m = 18 in Fig. |HJ We see - as in the case of the Kadanoff-Baym evolution - 
that all initial states show an approximately exponential decrease in time. The relaxation 
rates - as determined from the slope of the exponential decline - are also approximately 
the same for all distributions. They are given by T^ 1 « 0.98 ■ 1(T 2 for distribution Dl, 
T° 2 « 1.01 ■ 10~ 2 for distribution D2 and T^ 3 w 1.07 ■ 10~ 2 for distribution D3. Thus the 
relaxation rates are in the same range as those found within the full Kadanoff-Baym theory 
(Section IV. D in Ref. [HI])- This is exactly the result one expects from the analytical 
estimate for the chemical potential relaxation rate. In Section IV.D of Ref. [5T| we have 
found that the relaxation rate can be explained within a linearized evolution equation 
including only equilibrium properties, i.e. the equilibrium spectral and (Bose) distribution 
function. It is appropriate for small deviations from the equilibrium state in terms of the 
chemical potential /i. In the present case of the generalized transport equation (|39|) we 
encounter exactly the same situation. We know from the validity of the estimate that a 
linearized description is meaningful. Thus the evolution within the first order gradient 
equation should yield a comparable result as long as the equilibrium properties are about 
equal. This is indeed the case since the final temperatures for the various initial states are 
approximately the same (also compared to the Kadanoff-Baym case). They are given by 
T^/m « 1.819, T™/m « 1.827, T^/m w 1.834 compared to T eq /m « 1.835 in the full 
Kadanoff-Baym case. Consequently, the same similarity holds for the spectral function, 
which is determined in equilibrium by the temperature T and the coupling strength A. 
Therefore we can conclude, that the generalized transport equation (|39|) is sufficient to 
describe the correct relaxation of the chemical potential /x. 

Finally, we study the relaxation of the chemical potential as a function of the coupling 
strength A. To this aim we display in Fig. [§] the relaxation of the chemical potential 
\x for the initial distribution D2 for coupling constants A/m = 8, 10, 12, 14, 16, 18 as 
obtained form the generalized transport equation (|3*9*|) . For all coupling constants A the 
chemical potential /i is reduced exponentially in time oc exp(— ■ t) and thus allows for 
the determination of a proper relaxation rate T^. 

The results for the relaxation rates of the chemical potential as a function of the 
coupling strength A are displayed in Fig. [TUJ They are scaled by the coupling constant 
squared (A 2 ) in order to take into account the overall coupling dependence of the collisional 
self-energy. We see that the relaxation rate T M increases much stronger than quadratically 
with the coupling strength. Whereas the relaxation is very weak for small and medium 



21 



couplings X/m < 8, it increases considerably for larger interaction strength. By inspection 
of the analytical estimate for the relaxation rate of the chemical potential (Eq. (4.28) in 
Ref. [ST]) we find an explanation for this strong dependence: At first here an overall factor 
of the coupling constant squared (A 2 ) enters the expression for the relaxation rate, that 
stems from the collisional integral in terms of the scattering (sunset) self-energies. This 
factor, only, would yield a constant line in Fig. El and thus underestimate the observed 
behaviour significantly. Therefore, one has to keep in mind the additional A dependence of 
the (equilibrium) spectral function that strongly influences the estimate for the relaxation 
rate since it appears in the energy-momentum integration weights several times. Thus 
the relaxation rate is determined explicitly via the i^-term (Eq. (4.26) in Ref. (HIJ) and 
implicitly through the spectral functions in both contributions, K\ and K% (Eqs. (4.26) 
in Ref. j3TJ), by the coupling strength A in a nonlinear way. 

We conclude that - though there is a small relative shift of the different time scales of 
kinetic and chemical equilibration as a function of the coupling strength A with respect 
to the full KB solutions - the results of the generalized transport equations are very 
similar. The differences we attribute to higher order multi-particle effects in off-shell 
transitions. While the kinetic equilibration proceeds approximately with the coupling 
constant squared (as indicated by the calculations for non-polar-symmetric systems), the 
chemical relaxation rate is a higher order process in A as seen from Fig. Thus the 
chemical equilibration moves with increasing coupling strength A to earlier times relative 
to the kinetic relaxation, which is governed by terms ~ A 2 (see next Subsection). 

3.5 Quadrupole Relaxation 

In this Section we no longer restrict to polar symmetric systems and discuss the time 
evolution of more general initial distributions within the generalized transport approxi- 
mation (|39J). We start with conditions similar to those employed in Section V of Ref. [5T] . 
but combined with the initialization scheme for the semiclassical limit. Again - as in the 
full Kadanoff-Baym and the Boltzmann case (Section V of Ref. j3TJ) - the decrease of 
the quadrupole moment of the distribution 



is approximately exponential in time (oc exp(— Tq ■ t)) and thus allows for the extraction 
of a quadrupole damping rate Tq. The scaled quadrupole damping rates - as obtained 
for the two initial distributions dl and d2 (cf. Section V of Ref. (SJJ) - are displayed 
in Fig. ^JJas a function of the coupling strength X/m. The calculations show that the 
quadrupole relaxation rates within the semiclassical approximation ()39|) for both initial 
distributions dl and d2 is well within in the range of the full Kadanoff-Baym and the 
on-shell Boltzmann case. Additionally, the quadrupole relaxation rate is rather flat in the 
coupling A when divided by the coupling constant squared (A 2 ) as already observed for 
the other two evolution schemes in Section V of Ref. [HI]. Nevertheless, the relaxation 
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in the full KB-theory Q proceeds slightly faster than in the transport limits for large 
couplings. The latter effect is again attributed to higher order off-shell transition effects 
which are no longer incorporated in the generalized transport equation. 



3.6 Validity of the Gradient Approximation 

As we have seen in the previous Subsections the generalized transport equation (fHHJ) 
leads to a good agreement with the Kadanoff-Baym dynamics. This indicates that the 
semiclassical limit can be applied without loosing essential features of the full quantum 
dynamics for homogeneous systems (in particular for the initializations Dl, D2 and D3). 
We recall, that the underlying assumption for the validity of the first order gradient 
expansion scheme - which has been used to derive the generalized transport equation - is 
that all functions are slowly evolving in the mean space and time coordinates. Thus, in 
comparison to the first order time derivatives, the second order time derivatives should 
be small, such that they can be neglected to a good approximation. In this Subsection 
we will study now this criterion in a more quantitative way. 

To this aim we consider as a relative measure the energy-momentum integrals over 
the absolute value of first and second order time derivatives of various functions entering 
the generalized transport equation. Explicitly this measure is given at time t by 



for an arbitrary function F in Wigner space. In the following we take into account the 
time derivatives of the Green functions and the self-energies, i.e. 



The time evolution of these measures is shown in Fig. ^] for the first and the second 
order time derivatives of the Green function and in Fig. EH for the collisional and the 
retarded self-energies. The calculation has been performed for the initial distribution D2 
with a coupling strength of X/m — 18. For the Green functions as well as the self-energies 
the second order time derivatives are about one order of magnitude smaller than the first 
order expressions. Thus the underlying assumption of the first order gradient expansion 
is fulfilled very well indicating that the results obtained within the semiclassical scheme 
should match with those for the full quantum evolution to a large extent. This is exactly 
what we have found from the explicit comparison of the evolution of the equal-time Green 
functions as well as the spectral functions in Sections 13.21 and 13.31 

3.7 Generalized Transport in Botermans-Malfliet Form 

In this Subection we will perform a comparison of the generalized transport equation in 
the original Kadanoff-Baym (KB) form (|39|) with the modified Botermans-Malfliet (BM) 
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F G { iG < 1 iS < , Re G R , i?e E R } . 
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form (|45|). As discussed in detail in Section 2.2 the latter form results from the replacement 
of the collisional self-energy by 

iS< -> 2<5< -f/A (52) 

in the second Poisson bracket on the l.h.s. of the original kinetic equation (j39jl . This 
replacement leads to a consistent first order equation in the gradients and achieves con- 
sistency of the resulting transport equation with the corresponding generalized mass-shell 
relation (EEUl) . 

In Fig. E]we compare the time evolution within the generalized transport equation 
in the KB form to the consistent equation in BM form. In this respect several momentum 
modes of the equal-time Green function are displayed evolving in time from an initial 
distribution D2 for a coupling constant X/m — 16. We find that the deviations between 
both approximations (KB and BM) are rather moderate. Only for very small momentum 
modes \ p\/m < 1.6 deviations between both calculational modes are visible. For the very 
low momentum modes the range of difference starts at t • m ~ 10 and extends to t ■ m ~ 
100 for the non-zero modes. For the zero momentum mode the deviation lasts even longer. 
In this region the semiclassical transport in the BM form is slightly 'slower' than in the 
original KB choice. Nevertheless, also the BM form exhibits the typical overshooting 
behaviour of the low momentum modes beyond the stationary limit as observed for the 
KB form. However, the maxima are shifted slightly to later times. Finally, both gradient 
approximations converge in the long-time limit to very similar configurations. 

As the final part of the comparison we investigate the approximation (p)2*|) quantita- 
tively and introduce a measure for the absolute strength of the collision term by inte- 
grating the collision rates over energy and momentum. Since the collision term vanishes, 
if the substitution (|5^|) holds exactly, i.e. in equilibrium, the absolute size gives an idea 
about the validity of this replacement at the zero order level. Explicitly the measure reads: 



M col {t) 



d D p 1 
(2^p |2^T 



zS< iG > - zE> iG< 



(53) 



Furthermore, we define a measure for the deviation of the second Poisson bracket in the 
original and in the consistent Botermans-Malfliet formulation. It is given by the integra- 
tion over the absolute differences of the Poisson terms as 



M of f(t) 



d u p 1 



{iJ: < ,ReG R } - {iG < j, ReG R } 



(54) 



In Fig. E3we display the time evolution of the measure (|53|) (solid line) and the measure 
(I54J) (dashed line). The calculation has been performed for the initial state D2 with a 
coupling constant X/m = 16. We find that both measures decrease as a function of time 
in accordance with the equilibration of the system. Additionally, the contribution from 
the difference of the Poisson brackets (|54|) is always smaller than the one stemming from 
the collision term (|53|). This indicates that the replacement of the collisional self-energy 
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zS < — > iG < - T/A is more reliable when it takes place at the gradient level in accordance 
with the assumption of Botermans and Malfliet. However, the relative suppression is not 
very large. 

In summary we point out that the approximation of the full Kadanoff-Baym dynamics 
by the generalized transport equations in Kadanoff-Baym ()39|) or Botermans-Malfliet form 
f!45|) holds very well for the different momentum modes of the Green function iG < itself. 
Slight deviations are only visible for the zero momentum mode at early to intermediate 
times (Figs. E3 and IT^|) for a logarithmic representation of the time axis. Consequently, 
the characteristic features of quantum equilibration obtained for the full Kadanoff-Baym 
theory are retained in the generalized transport limits. The validity of these transport 
equations - based on a first order gradient expansion - could be shown explicitly, since 
second order gradient terms turned out to be smaller by more than an order of magnitude 
(cf. Figs. dandHSJ). 

4 Summary and Outlook 

In this work we have studied the quantum time evolution of 4 -field theory for homo- 
geneous systems in 2+1 space-time dimensions for far-from-equilibrium initial conditions 
in extension of our studies in Ref. jSI]. The three-loop approximation for the CTP 2PI 
effective action has been employed, i.e. the tadpole and sunset self-energies. The tad- 
pole contribution corresponds to a dynamical mass term whereas the sunset self-energy 
is responsible for dissipation and an equilibration of the system. Since both self-energies 
are ultraviolet divergent, they had to be renormalized by including proper counter terms. 
The numerical solutions for different initial configurations out of equilibrium (with the 
same energy density) show, that the asymptotic state achieved for t — > oo is the same for 
all initial conditions in the full Kadanoff-Baym theory. 

In order to improve the standard on-shell (or quasiparticle) Boltzmann transport the- 
ory we have derived generalized off-shell transport equations from the Kadanoff-Baym 
equations in phase-space representation by restricting to first order derivatives in x and 
p. In fact, it could been shown explicitly, that second order derivatives (of self-energies 
and Green functions) are smaller than the first order derivatives by at least an order of 
magnitude. As a consequence the dynamics within the generalized transport formulation 
shows a very similar structure as the full quantum solution. This is clearly seen for the 
propagation of the equal-time momentum modes of the Green functions and the evolution 
of the spectral function for all configurations considered. It includes, in particular, the 
overshooting behaviour of the low momentum equal-time modes, which occurs at interme- 
diate times depending on the initial distribution. Furthermore, the evolution within the 
generalized transport equation leads to a stationary state in the long time limit, which 
exhibits a full off-shell equilibration with vanishing chemical potential. Even the relax- 
ation rates of the chemical potential obtained from the semiclassical evolution agree very 
well with those of the full Kadanoff-Baym theory. The dependence of the chemical relax- 
ation rate on the coupling constant is rather non-trivial, since it is strongly affected by 
the equilibrium spectral function. Thus we conclude that the inclusion of the dynamical 
spectral function - as inherent in the semiclassical approximation of the KB equations - 
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surpasses the shortcomings of the on-shell Boltzmann limit. 

Moreover, we have shown that the generalized transport equation in Kadanoff-Baym 
(KB) form (|39|) and in Botermans-Malfliet (BM) form ()45|) lead to comparable results for 
the time evolution of the initial configurations considered. Only for the time-dependent 
occupation of low momentum modes slight differences have been observed. This is a typi- 
cal quantum phenomenon related to the large de Broglie wavelength of the low momentum 
modes. We stress that the generalized transport equation in BM form can be used for the 
off-shell description of realistic heavy-ion collisions El HH IS] - as well as many other 
interacting systems in various fields of physics - since it allows for an economic solution 
within an extended testparticle ansatz. 

Coming back to the questions raised in the Introduction, we have to answer the first 
one by 'probably no' since off-shell quantum transitions in comparison to on-shell two- 
body scattering only lead to slightly smaller relaxation times e.g. for the quadrupole 
moment in momentum space. On the other hand, the off-shell transitions are important 
(in case of the 4 -theory) for chemical equilibration, which essentially proceeds via 1 «-> 3 
transitions jHIj. The latter are forbidden in the on-shell quasiparticle limit due to energy 
and momentum conservation. We recall that also 0^4 processes are incorporated in 
the full Kadanoff-Baym theory though their effect was found to be very limited for the 
dynamical configurations investigated here. 

Nevertheless, this still leaves us with the open question on the relative importance 
of many-body scattering processes n m for n > 2, m > 2. Here, an appropriate 
formulation of n <-> m off-shell and on-shell scattering processes has already been given 
in Ref. [SH] as well as an application to the problem of antibaryon annihilation and 
recreation by several mesons in relativistic nucleus-nucleus collisions. A related study 
on the problem of 2 <-> 3 transitions (e.g. bremsstrahlung for massless particles) on the 
parton level will be presented in the near future [60 . 
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Figure 1: Initial momentum distributions for the initializations Dl, D2 and D3. For the 
polar symmetrical settings in momentum space the distributions are shown as function of 
p x (for p y = 0). 
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Figure 2: Initial equal-time Green functions for the initializations Dl, D2 and D3. The 
equal-time Green functions are displayed for polar symmetric systems as a function of p x 
(for p y = 0) with coupling constant X/m = 18. 



30 



Q_ 

v 

CD 



0,4 
0,2 
0,0 



1 1 1 — I — I I I I I 1 1 1 — I — I I I I 

KB-transport 



2,2 
2,0 
1,8 
1,6 
1,4 
1,2 
1,0 
0,8 

o,6 ^ ; 



- 1 — i — i — i — 1 1 1 1 



— - distribution D1 

distribution D2 

distribution D3 



Urn = 18 

p/m = O.CT 



"■■■■■■»■» s ■ = , ~ JXi — . 



p/m =.0.8; 

p/m = 1 .6^ 
p/m = 2.4^ 

~p/m~= 3]£ 



- 1 1 1 — 1 — 1 — 1 1 1 J 1 1 1 — 1 — 1 — 1 1 1 J 1 1 1 — 1 — 1 — 1 1 1 J 

10 100 1000 



time t m 



Q_ 

v 

CD 



2,2 
2,0 
1,8 
1,6 

1,4- 

1,2- 

1,0- 

0,8 

0,6 

0,4 

0,2 



0,0- 



1 1 1 — 1 — 1 1 1 1 | 1 1 1 — 1 — 1 1 1 1 | 1 1 111 rrrj 



KB-theory 



distribution D i 

distribution D2 

distribution D3 



Xlm = 18 

p/m=0.0- 

p/m=0.8- 



J-..- >^.^a-' 



p/m=1.6 _ 
_p/m=2jr 

p/m=3.2- 



-I 1 1 — 1 — 1 1 1 1 | 1 1 1 — 1 — 1 1 1 1 | 1 1 1 — 1 — 1 — 1 1 1 | 



1 



10 



100 



1000 



time t m 



Figure 3: Upper part: Evolution of several momentum modes \ p\/m = 0.0, 0.8, 1.6, 
2.4, 3.2, 4.0 of the equal-time Green function on a logarithmic time scale for the different 
initializations Dl, D2 and D3 for the generalized transport equation (J39)) . Lower part: 
same as above but for the full Kadanoff-Baym equation Q . 
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Figure 4: Time evolution of the spectral function A(p, po,t) for the initial distributions 
Dl, D2 and D3 (from top to bottom) with coupling constant X/m = 18 for the two 
momentum modes \ p\/m = 0.0 (l.h.s.) and \ p\/m = 0.0 (r.h.s.). The spectral function 
from the transport eq. (JH is shown for times t ■ m = 5, 20, 60, 120, 240, 360, 480, 600 
as indicated by the different line types. 
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Figure 5: Time evolution of the on-shell energies u(p, t) of the momentum modes \p\/m = 
0.0 and \p\/m = 2.0 for the different initializations Dl, D2 and D3 with X/m = 18 in 
the semiclassical KB limit (|39|) (l.h.s.). The on-shell self-energies are extracted from the 
maxima of the time-dependent spectral functions. The respective results from the full 
Kadanof-Baym theory are displayed on the r.h.s. 
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Figure 6: Time evolution of the on-shell widths — JmE R (p, ou(j>, t),t)/u(p, t) of the mo- 
mentum modes \ p\/m = 0.0 and \ p\/m = 2.0 for the different initializations Dl, D2 and 
D3 with A/m = 18 in the semiclassical transport eq. (|39|) (l.h.s.). The respective results 
from the full Kadanof-Baym theory are displayed on the r.h.s. 
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Figure 7: Spectral function A for various momentum modes as a function of energy 
Po/m at the late time t ■ m — 600 (lower part) for initial distribution D2 with coupling 
constant A/m = 18 in the semiclassical transport limit ()39|) (l.h.s.). The corresponding 
distribution function N (at the same time for the same momentum modes) is shown in the 
upper part. All momentum modes can be fitted for all energies by a single Bose function 
of temperature T/m = 1.827 and a chemical potential fi/m close to zero. The respective 
results from the full Kadanof-Baym theory are displayed on the r.h.s. 
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Figure 8: Relaxation of the chemical potential fi/m for the three initial distributions Dl, 
D2, D3 within the generalized transport equation (JBTJj) in a logarithmic repesentation as 
a function of time. The decrease is approximately exponential in time oc exp(— ■ t) and 
the relaxation rate is practically equal for all initializations. 
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Figure 9: Relaxation of the chemical potential /i/m for initial distribution D2 with 
coupling constants X/m = 8, 10, 12, 14, 16, 18 using the generalized transport equation 

dm. 
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Figure 10: Scaled relaxation rate T M /A 2 of the chemical potential as a function of the 
coupling strength X/m for the initial distribution D2. 
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Figure 11: Scaled Relaxation rate for the generalized transport equation as a function of 
the coupling strength (half-filled symbols) for the initial distributions dl (squares) and d2 
(triangles). Additionally the results obtained within the Kadanoff-Baym (full symbols) 
and the Boltzmann calculation (open symbols) are shown for comparison. 
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Figure 12: Time evolution of the first /second order time derivative measure for the Green 
functions iG K and ReG R for initial distribution D2 with coupling constant A/m = 18. 
During the whole evolution the second order time derivatives are more than an order of 
magnitude smaller than the first order derivatives supporting the validity of the gradient 
expansion. Note that time derivatives in our convention are derivatives with respect to 
the dimensionless quantity tm. 
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Figure 13: Time evolution of the first/second order time derivative measure for the 
functions zS < and ReY> R for initial distribution D2 with coupling constant A/m = 18. 
During the whole evolution the second order time derivatives are significantly smaller than 
the first order contributions. Note that time derivatives in our convention are derivatives 
with respect to the dimensionless quantity tm. 
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Figure 14: Time evolution of various equal-time Green function momentum modes 
within the generalized transport equation (original KB form, solid lines) and within the 
Botermans-Malfiet (BM) form (dashed lines) for initial distribution D2 with coupling 
constant X/m = 16. 
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Figure 15: Time evolution of the measure functions for the collision term (|5H|) and the 
difference of the second Poisson bracket in standard (KB) and consistent (BM) form (|H%j) 
for the initial distribution D2 with coupling constant A/m = 16. Both contributions 
decrease in time in accordance with the equilibration of the system. The Poisson contri- 
bution (dashed line) is always significantly smaller than that of the collision term (solid 
line) . 
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